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Abstract. The electromagnetic response of subwavelength wires coated with a 
graphene monolayer illuminated by a linearly polarized plane waves is investigated. 
The results show that the scattering and extintion cross-sections of the coated wire 
can be dramatically enhanced when the incident radiation resonantly excites localized 
surface plasmons. The enhancements occur for p-polarized incident waves and for 
excitation frequencies that correspond to complex poles in the coefficients of the 
multipole expansion for the scattered field. By dynamically tuning the chemical 
potential of graphene, the spectral position of the enhancements can be chosen over a 
wide range. 


PACS numbers: 78.67.Wj, 73.20.Mf, 78.20.Ci, 81.05.ue, 73.21.-b, 78.68.+m 

1. Introduction 

Advanced light manipulation relies heavily upon controlling the fundamental processes 
of light-matter interaction. This paper is concerned with light scattering, a process 
which plays a key role in current light control applications like light trapping [Tj , 
Anderson localization 12113] or photonic band gaps [3]. Particularly, light scattering 
by subwavelcngth particles has received considerable attention in relation to plasmonic 
systems, where some exciting experimental and theoretical developments can be 
produced when incident radiation is coupled to the collective electron charge oscillations 
known as localized surface plasmons j5ll6li71l8ll9] . Excitation of localized surface plasmons 
by an incident electromagnetic held at the frequency where resonance occurs results 
in an enhancement of the local electromagnetic held as well as in the appearance of 
intense absorption bands. Moreover, the efficiency of light scattering, determined by the 
strength of the electric held in the material, is also enhanced under resonant excitation 
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of surface plasmons. The fact that the details of the resonance are highly sensitive to 
particle characteristics such as material, size, shape and enviroment constitutes the 
precise property which has prompted an exponential growth in the exploitation of 
localized surface plasmons in nanoscale optics and photonics [7], giving applications 
such as molecular-recognition elements [8], amplifiers [9] and biosensors. Applications 
and new developments in the field of light scattering by plasmonic particles have been 
recently reviewed in [5]. 

Natural plasmonic materials such as metals or metal-like materials are nonmagnetic 
(magnetic permeability of free-space) and have electric permittivity with a negative 
real part. Currently available plasmonic materials have a number of shortcomings, 
most notably power losses and bad tunability, which limit the further development of 
plasmonics. The recent advent of graphene -a monolayer of carbon atoms arranged 
in a hexagonal lattice- have just meet the need of surface plasmons since it offers 
relatively low loss, high confinement and good tunability, three characteristics that make 
graphene a promising plasmonic material alternative to noble metals [TOlIl 111121113] . The 
waveguiding characteristics of graphene plasmons have been investigated in different 
planar structures, such as infinite sheets ra, ribbons PS 03] or grooves H3- Since 
the wave vector needed to resonantly excite these plasmons is always greater than that 
available to incident electromagnetic radiation, the electromagnetic response of graphene 
surface plasmons has been investigated for phase matching configurations such as the 
prism [18] and grating [191120] I2T] couplers as well as for layered media [22,23]. 

Quite recently non-planar configurations [211125112611271128] have also begun to attract 
particular attention since cylindrical and spherical structures have been experimentally 
constructed for a number of applications [29], !30, :3T[,[32], [33[[ 34] . The electromagnetic 
response and the plasmonic properties of graphene-coated nanospheres have been 
recently examined in [28] and the modal characteristics of surface plasmons propagating 
along the axis of graphene-coated wires have been studied in [2312511261127] . However, the 
resonant excitation of localized plasmon modes in graphene-coated wires has not been 
investigated yet. The primary motivation of the work described here is to study the 
plasmonic response of graphene-coated 2D particles (wires). For metallic or metallic- 
like substrates, the coating is expected to modify the localized surface plasmons already 
existing in the bare particle. On the other hand, for substrates made of nonplasmonic 
materials the coating is expected to introduce localized surface plasmons mechanisms 
which were absent in the bare particle. In both cases, the good tunability of the 
graphene coating is expected to lead to unprecedented control over the location and 
magnitude of the particle plasmonic resonances. In this paper we study the scattering 
of electromagnetic waves normally incident on graphene-coated cylinders with circular 
cross-section, and present results for both dielectric and metallic substrates. 

The plan of the paper is as follows. In Section [2] we sketch an analytical method of 
scattering based on the separation of variables approach and obtain a solution for the 
scattered fields in the form of an infinite series of cylindrical multipole partial waves. 
In Section [3] we give examples of scattering and extinction cross-section corresponding 
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to subwavelength wires tightly coated with a graphene layer. Both transparent and 
metallic susbstrates are considered. In both cases, great enhancements are observed for 
p-polarized incident radiation and for excitation frequencies corresponding to localized 
plasmons of the graphene-coated cylinder. Finally, in Section [4] we summarize and 
discuss the results obtained. The Gaussian system of units is used and an exp (—i cot) 
time-dependence is implicit throughout the paper, with to the angular frequency, t the 
time, and i = \/—T. 

2. Theory 


/ 



Figure 1 . Schematic of the scattering problem 

We consider a graphene-coated cylinder with circular cross-section (radius R) 
centered at x=0, y —0 (see Figured]). The wire substrate may be dielectric or conducting 
(electric permittivity and magnetic permeability /ii) and the coated wire is embedded 
in a transparent medium (electric permittivity £2 and magnetic permeability // 2 ). The 
graphene layer is treated as an inhnitesimaly thin, local and isotropic two-sided layer 
with surface conductivity cr((u) given by the Kubo formula [35][36j. The wave vector of 
the incident radiation is directed along x. In this case, the scattering problem can be 
handled in a scalar way since the solution to any incident polarization can be described 
as a linear combination of the solutions obtained in two fundamental scalar problems: 
electric field parallel to the main section of the cylindrical surface (p polarization or 
E ~=0 modes) and magnetic field parallel to the main section of the cylindrical surface 
(s polarization or H z =0 modes). 

To obtain analytical solutions to the scattering problem we closely follow the usual 
separation of variables approach [37jl38]. In a first step, the non-zero components of 
the total electromagnetic field along the axis of the cylinder for each polarization case, 
denoted by F(p, tp), are expanded as series of cylindrical harmonics, one for the internal 
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region (p < R , subscript 1) and another one for the external region (p > R , subscript 
2). When the incident electric held is contained in the x — y plane (p polarization), the 
expansions for the total magnetic held along the axis of the cylinder are 

OO 

Fi(p,<p)= c n J n{k\p) exp imp, (1) 

n =—oo 

OO 

F 2 (p, ( p)= ^2 [AoFJnihp) + a n H^(k 2 p)] exp imp, ( 2 ) 

n =—oo 

where a n and c n are unknown complex coefficients, kj = p^/sjp] (j = 1,2), c is the 
speed of light in vacuum, A 0 is the amplitude of the incident magnetic held (parallel to 
z) and J n and > are the n-th Bessel and Hankel functions of the hrst kind respectively. 
In the same manner, when the incident magnetic held is contained in the x — y plane 
(s polarization), the expansions for the total electric held along the axis of the cylinder 
are 

OO 

F\ (p, <p) = ^2 d nJn(kip) exp imp,, (3) 

n=— oo 
oo 

F 2 (p, <p) = ^2 [Bo i n Jn(k 2 p) + b n H^\k 2 p)] expimp , ( 4 ) 

n=— oo 


where B 0 is the amplitude of the incident electric held (parallel to z) and b n and d n 
unknown complex coefficients. 

In a second step, the boundary conditions at p = R are invoked to obtain the 
unknown coefficients a n , b n , c n and d n in terms of the incident amplitudes A 0 and 
Bo. The boundary conditions for the graphene-coated cylinder are: i) the tangential 
component of the total electric held E is continuous and ii) the discontinuity of the 
tangential component of the total magnetic held H is proportional in magnitude to the 
surface current density (that is, to the graphene surface conductivity). In terms of F, 
the boundary conditions at p = R can be expressed as 


1 d F 2 _ 
£ 2 dp 

1 dF\ 

£ i dp 

Am dF 1 

and F 2 — F x =-cx—— , 

U£]_ Op 

(5) 

for p polarization, and 




Fi = F 2 

and 

1 3F 2 1 d Fi Amu 

/i -2 dp pi dp c 2 1 

(6) 


for s polarization. Finally, the amplitudes of the scattered held can be written as 


-i n £xk 2 Jn{xx)J' n {x 2 ) - + ^-a^ikifaJ^x^J^xt) 

£ik 2 J n (xi)H' n {1 \x 2 ) - £ 2 kiJ' n {x l )Hn\x 2 ) + ^-a^ikxfa J' n {xi)H'^ l) {x 2 ) 


( 7 ) 


~i n £ 2 k 1 J n (x 1 )J , n (x 2 ) - £ik 2 J' n {x 1 )J n (x 2 ) + ^-a^-ik 1 k 2 J n (x 1 )J n (x 2 ) 

_ ; __ —B 

£ 2 kiJ n (x 1 )H , ^ 1 \x 2 ) - £ik 2 Jn(xi)Hn\x 2 ) + ^ikik 2 J n {xi)Hn\x 2 ) 
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while the amplitudes of the fields inside the wire are given by 


£ik 2 i n Jn{x 2 )H^ l \x 2 ) - J' n {x 2 )H^\x 2 ) 



e 2 kii n J n {x 2 )Hn l \x 2 ) - J' n (x 2 )Hn\x 2 ) 



where X\ = k\R, x 2 = k 2 R. 

As in the problem of scattering by a wire without a graphene coating, the amplitudes 
given by equations ©-(HOD allow us to obtain the electromagnetic field everywhere 
in space as well as other quantities of interest such as differential, absorption and 
extinction cross-sections (see [0711381 for details). The multipole coefficients a n and 
b n for the scattered field have essentially the same form as those corresponding to a 
bare wire [37[l38[. except for additive corrections proportional to a in numerator and 
denominator. Similar additive corrections have been recently reported for graphene- 
coated spheres [28]. 

3. Resonant excitation of localized surface plasmons 

To simulate the electromagnetic response of graphene-coated cylinders we assume that 
cr(cn), the complex surface conductivity of graphene, is given by the high frequency 
expression derived from Kubo formula (equation (1), Ref. [35]). including interband and 
intraband transition contributions. Apart from the angular frequency ui, the value of a 
depends on the chemical potential p c (controlled with the help of a gate voltage), the 
ambient temperature T and the carriers scattering rate 7 C . For the intraband (cr mtm ) and 
interband (cr mter ) contributions to ct(oj) we have used the following expressions [35,36] 



( 11 ) 


(cu + 2/i c ) 2 


(u ~ 2/i c ) 2 + (2T) 2 


where h is the reduced Planck constant and e the elementary charge. In all the 
calculations we have used T = 300° K and 7 C = 0.1 meV. 

In Figure [2] we plot the scattering and absorption cross-sections (per unit length) 
for a wire with a radius R = 0.5 /mi, made with a nonplasmonic, transparent material 
(£1 = 3.9, pi = 1) in a vaccum (p 2 = £2 = 1)- The excitation frequencies are in the 
range between 5 THz (incident wavelength 60 pm) and 30 THz (incident wavelength 
10 /mi), the incident wave is p-polarized and different values of p c have been considered. 
In Figure [2^ the scattering cross-section spectrum of the bare wire (continuous line) 
is also given as a reference. The absorption cross-section of the bare wire vanishes 
identically and is not shown in Figure [2}o. We see that while the scattering cross-section 
spectrum of the bare wire does not show any plasmonic feature in this spectral range, 
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Figure 2. Scattering (a) and absorption (b) cross-section spectra for a graphene- 
coated cylinder illuminated by a p-polarized plane wave and for different values of 
the chemical potential g c . The scattering cross-section curve corresponding to the 
uncoated cylinder is given as a reference. The incident wave is p-polarized, R = 0.5 /im, 
£\ = 3.9 and gi = ji 2 = £2 = 1. 
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great enhancements in both the scattering and the absorption cross-section spectra 
occur at wavelengths near 46.84 /irn (/z c = 0.3 eV), 36.35 /irn (ji c = 0.5 eV) and 
28.81 /irn (/i c = 0.8 eV). At these spectral positions we obtain enhancement factors in 
the scattering cross-section of around four orders of magnitude. Other local maxima, 
particularly noticeable in the absorption cross-section spectra, occur at wavelengths near 
46.84 /irn, 33.06 /irn (for /x c = 0.3 eV), 36.35 /irn, 25.64 /irn, 20.92 /irn (for /x c = 0.5 eV) 
and 28.81 /irn, 20.29 /irn, 16.55 /irn (for /i c = 0.8 eV). The local maxima observed in the 
scattering and absorption cross-sections displayed in Figure [2] correspond to plasmonic 
resonances of the graphene-coated wire and can be obtained from the complex poles of 
the coefficients a n and c n of the multipole expansions for the electromagnetic field. The 
correspondence can be clearly seen by noting that the spectral positions of the peaks are 
in very good agreement with the real part of the complex wavelengths A n given in Table 
1 (n = 1, ...4), with A n representing the complex root of the common denominator 
g n { A) in equations (J7J) and (j9]) , and A = 27 tc/uj. For the numerical calculation of A„ 
we have used a Newton-Raphson method and the same parameters used for curves in 
Figure [2j We observe that the scattering cross-sections in Figure [2^ display minima 
at wavelengths near 36.09 /irn (/z c = 0.3 eV), 28.04 /irn (/r c = 0.5 eV) and 22.27 /irn 
(/i c = 0.8 eV). We found that the spectral position of these minima are in excellent 
agreement with the real part of the complex wavelengths where the numerator of the 
coefficient cq becomes zero. Taking into account that a minimization of the scattering 
cross-section is a necessary condition to obtain invisibility or transparency of an object, 
and that the scattering cross-section of the bare wire does not exhibit a minimum 
in this spectral range, the observation of these minima is relevant in the context of 
invisibility cloacks. These results, together with those shown in Figure [2l illustrate a 
great advantage of the graphene coating, which provides unprecedented control over the 
spectral location of the resonances of the wire via changes in the chemical potential /i c , 
ultimately controlled through electrostatic gating. 

The correspondence between plasmonic resonances of the graphene coating and 
the enhancement of scattering and absorption in the coated wire is further evidenced 
when the cross-sections of the same systems considered in Figure [2] are calculated 
for s-polarized -instead of p-polarized- illumination. In this case (not shown), no 
enhancements in the cross-sections of the graphene-coated wire are observed, which 
is consistent with the fact that localized surface plasmons are not supported in this 
polarization mode, since the electric held in the graphene coating can only induce 
electric currents directed along the wire axis, and not along the azimuthal direction 
ip, a necessary condition for the existence of localized surface plasmons in the graphene 
circular cylinder. 

In Figure |3k we plot the spatial distribution of the electric held normalized to the 
incident amplitude for the wire considered in Figure [2] and p c = 0.8 eV. The incident 
wavelength is A = 28.809 /irn, a value equal to the real part of the complex pole Ai (see 
Table 1) of the multipole coefficients a\ and C\ and for which the strongest maxima in 
the scattering and absorption cross-sections for /i c = 0.8 eV occur. We observe that the 
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Figure 3. Map of | E(p, ip )| for a graphene-coated cylinder illuminated by a p-polarized 
plane wave. The incident wavelength is A = 28.809 /im (top), A = 20.295 /im (center) 
and A = 16.553 /xm (bottom). The chemical potential is /x c = 0.8 eV and all other 
parameters are the same as in Figure [2] The circle (continuous line) indicates the 
graphene layer. The top inset illustrates schematically the angular distribution of 
scattered power, the medium inset illustrates the values of |c n |, n = 0, 1, 2, 3 (the 
multipole coefficients inside the cilynder) and the bottom inset illustrates the values 
of |a n |), n = 0, 1, 2, 3 (the multipole coefficients outside the cilynder). 
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n 

p c = 0.3 eV 

p c = 0.5 eV 

p c = 0.8 eV 

1 

2 

3 

4 

46.8378 + *0.1216 

33.0651 + f0.0439 

26.9865 + z0.0292 

23.3660 + *0.0219 

36.3464+ f0.0961 
25.6361+ f0.0264 

20.9175 + *0.0175 

18.1086 -H0.0131 

28.8092 + f0.0875 
20.2956 + *0.0167 

16.5531+ f0.0109 

14.3272 -H0.0082 


Table 1: Complex roots A n of g n ( A) = 0 in , where g n is the common denominator of 
the coefficients a n and c n and A = 2nc/co. The parameters correspond to those used to 
calculate the scattering and absorption cross-sections shown in Figure [2j 


picture is very similar to that obtained for subwavelength metallic wires: the near fields 
are clearly of an electric dipole nature, with the dipolar moment oriented in the direction 
of the incident electric held (y). However, contrary to the case of subwavelength metallic 
(impenetrable) wires, where the held cannot penetrate into the interior regions and is 
limited to a surface layer of approximately one skin depth thick, in the graphene-coated 
dielectric wire there are regions where the interior held can be greatly enhanced. Taking 
into account that the electric held induced in the scatterer is much stronger than that 
of incident radiation (see scale in Figure Eb), the far-held intensities are dominated by 
an electric dipole pattern, as indicated in the inset in Figure |3h, where we have sketched 
the angular distribution of scattered power for this excitation wavelength. When the 
incident wavelength is A = 20.295 /mi, the scattering cross-section for p c = 0.8 eV (see 
Figure [2]) exhibits a very weak peak while the absorption cross-section for the same 
value of p c exhibits a strong enhancement (an enhancement factor near three orders of 
magnitude greater than the non-resonant case). At this wavelength, the electric held 
induced in the scatterer, although not as strong, is again stronger than that of incident 
radiation and of a quadrupolar, not dipolar, nature. This is shown in Figure Eb, where 
we plot the map of | E(p, ip) | normalized to the incident amplitude and calculated for the 
same parameters as in Figured, except the incident wavelength A = 20.295 /tm. Since 
the electric held induced in the scatterer is much stronger than that of incident radiation, 
the far-held intensities are dominated by an electric quadrupole pattern, as indicated 
in the inset in Figure Eb- In Figure Eh we show the map of \E(p,cp)\, normalized to 
the incident amplitude, for the same parameters as in Figure Eh and Eb, except that 
now the incident wavelength is changed to A = 16.553 /tm. For this value of A and for 
p c = 0.8 eV a weak peak can be observed in the absorption cross-section, although no 
peak is observed in the scattering cross-section. We observe that the spatial distribution 
of the electric held induced in the scatterer at this excitation wavelength is hexapolar- 
shaped. However, the far-held intensities are dominated by an electric dipole pattern, 
as indicated in the bottom schematic diagram in Figure Eb) which illustrates the values 
of \a n \ and \c n \, n — 0,1, 2, 3. Similar results have been obtained for the other values 
of chemical potential considered in Figure El 
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Figure 4. Extinction cross-section in a graphene-coated dielectric wire in vacuum 
{fi 2 = £2 = 1), explored as a function of the wire radius R and the excitation wavelength 
A. The dashed lines indicate the resonance position, numerically obtained as the poles 
of the coefficients a n . The incident wave is p-polarized, g c = 0.8 eV, £\ = 3.9 and 
Pi = 1 . 


To investigate the size dispersion of the localized plasmonic modes of graphene- 
coated wires, we show in Figure []] a color map of the extinction cross-section in the 
R — A plane. The incident wave is p-polarized, fi c = 0.8 eV, the ambient medium is 
vacuum, £4 = 3.9 and /p = 1. We observe that the enhancements in the extinction 
cross-section follow the dashed lines labeled with n. These lines indicate the resonance 
position, numerically obtained as the real part of the pole of the coefficient a n given by 
equation (17]). 

The examples so far highlight the attractive plasmonic features that a graphene 
coating can produce in intrinsically nonplasmonic wires. Next, we consider graphene- 
coated wires made of metallic or metallic-like substrates, that is, intrinsically plasmonic 
wires, where the localized surface plasmons of the graphene coating are expected to 
modify the localized plasmons already existing in the bare particle. We assume that the 
dispersive behavior of the interior electric permittivity is described by the Drude model 

,2 


£i(u) = £oo - 




UJ 2 + 


(13) 


with Sqq the residual high-frequency response of the material, u p the metallic plasma 
frequency and 7 m the optical loss-rate of the Drude material. This model is applied 
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Figure 5. Scattering cross-section spectrum in vacuum for a graphene-coated metallic 
wire with a radius R = 50 nm illuminated under p polarization. <r(w) is given by the 
Kubo expression with T = 300° K, y c = 0.1 meV and p c = 0.5eV. £i(w) is given 
by equation (fl3l) . with £<*, = 1, = 0.01 eV and Hui p = 0.5eV, 0.7 and 0.9eV. The 

scattering cross-sections for the bare metallic cylinders and for an empty graphene 
cylinder in vacuum (black curve, £i = £2 = 1 and /ii = /i 2 = 1) are also given as a 
reference. 


to describe strongly doped semiconductors [39] which allow dynamic manipulation of 
carrier densities and have plasma frequencies in the same range that the realizable Fermi 
energies of graphene. Here the situation is similar to that reported for graphene-coated 
Drude spheres [28] in the sense that the region of induced charges in the graphene 
coating is not separated from the region of induced charges in the metallic substrate. 
For the case considered in this paper, both regions coincide with the cylindrical surface 
of radius R and therefore a single hybridized plasmon is expected, rather than two 
(bonding and antibonding) hybridized plasmons with different energies, characteristic 
of systems where induced charges occur in two separated regions. In Figure |5] we plot the 
scattering cross-section spectrum in vacuum for a graphene-coated metallic wire with a 
radius R = 50 nm illuminated under p polarization. The interior electric permittivity 
£1 (oj) is given by equation (TT3]l with e 00 = 1, 7 m — 0-01 eV and different values of frw p 
(0.5 eV, 0.7 eV and 0.9 eV), while the values of the complex surface conductivity of 
graphene cr(co) have been obtained from the Kubo expression with T = 300° K, 7 C = 0.1 
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meV and /x c = 0.5eV. The scattering cross-sections for the bare metallic cylinders and 
for an empty graphene cylinder in vacuum (black curve, £i — e 2 — 1 and pi = fi 2 = 1) 
are also given as a reference. The excitation frequencies are in the range between 23.87 
THz (incident wavelength 12.57 urn) and 191 THz (incident wavelength 1.57 /ini). 

The curves in Figure [5] show that the hybridized resonances are always blueshifted 
compared to the resonances of the bare wire. This is consistent with the fact that 
the net effect of the graphene coating is to increase the induced charge density on the 
cylindrical surface of the metallic wire. The hybridized resonances can be controlled 
with the help of a gate voltage, which ultimately determines the value of the chemical 
potential p c . To explore the tunability of the hybridized resonances, in Figure |6] we 
give scattering cross-section spectra for a graphene-coated metallic wire ( hw p = 0.5eV) 
illuminated under p polarization for different values of the chemical potential /i c (0.3eV, 
0.5, 0.8 and 1.1 eV). 
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Figure 6. Scattering cross-section spectrum in vacuum for a graphene-coated metallic 
wire (huip = 0.7eV) illuminated under p polarization for different values of the chemical 
potential p c (0.3eV, 0.5, 0.8 and 1.1 eV). 


4. Summary 

Fueled by recent experiments [29ll30ll3Tll32] showing that, thanks to the van der Waals 
force, a graphene sheet can be tightly coated on a fiber surface, in this paper we have 
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investigated the electromagnetic response of circular cross-section wires coated with a 
graphene monolayer. The emphasis has been centered around the plasmonic properties 
of the coated wire, a system that not only is interesting for several applications, but 
provides a simple canonical model for the plasmonic characteristics of other graphene- 
coated 2D particles. We developed an analytical method of scattering based on the 
separation of variables approach and obtained a rigorous solution in the form of an 
infinite series of cylindrical multipole partial waves. We presented examples for both 
dielectric (nonplasmonic) and metallic (plasmonic) wire substrates. In the first case the 
graphene coating introduces surface plasmons mechanisms which were absent in the bare 
particle. In the second case, there is an interaction between the plasmons supported 
independently by the bare metallic wire and those suported by the graphene coating. 
In both cases, the examples show that the good tunability of the graphene coating leads 
to unprecedented control over the location and magnitude of the particle plasmonic 
resonances. 
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